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EXTRACTING FINE-TUNED ESTIMATES FROM CORRELATION 
FUNCTIONS EVALUATED AT A LIMITED NUMBER OF VALUES 
FIELD OF THE INVENTION 

The present invention relates to signal processing and, more particularly, to tech- 
5 niques for extracting fine-tuned estimates from correlation functions evaluated at a 
limited number of values. 
BACKGROUND OF THE INVENTION 

In the context of signal processing, when a signal is transmitted to a receiving 
device, there is a certain amount of unknown delay associated with the propagation 

10 of the signal. The original signal is not received instantaneously at the receiving 
device. A calculated estimate of the delay associated with the received signal is 
generally limited to the resolution of the sampling rate of the received signal For 
example, if the received signal can be represented as a continuous integrable function, 
x(t), then the received signal, which typically contains noise, may be converted to a 

15 discrete sequence of values in time through sampling. The discrete sequence of values 
is herein referred to as sampled data. As the sampling rate increases, the sampled 
data becomes a more accurate representation of x(t). Due to the finite sampling rate, 
the delay value that is calculated from the sampled data is not as precise as the delay 
value calculated from a continuous function. 

20 In one approach, a more accurate delay value estimate may be obtained by ad- 
justing the sampling times. For example, correlation calculations can be performed 
for a pair of delay values that produce equal values in the magnitude calculations of 
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the correlation integrals. If the magnitude values corresponding to the pair of delay- 
values are high enough, then it may be assumed that the peak magnitude, and hence 
the true time delay, lies halfway between the pair of delay values (used in the delay 
lock loop). However, a significant drawback with this approach is that it assumes 
5 that the sampling times can be adjusted, and that it is already known where to look 
approximately. If the sampling times cannot be adjusted, (i.e., the set of sampled 
data with which to work is fixed), then in general a different reference signal needs 
to be used in computing the correlation integrals at the different delay. It may not 
be enough to just time shift the reference by an integer number of samples. 

10 Another approach for estimating a delay value at a level of precision that is higher 
than the precision corresponding to the sampling rate, is to adopt a fine-grid for 
a hypothesized range of delay values. In such an approach, the sampled data is 
correlated against the reference signal by shifting the reference signal according to 
the fine gridding of delay values. The correlation between the received signal and 

15 the reference signal is performed by calculating the In Phase ("I") and Quadrature 
("Q") correlation integrals for the range of hypothesized delay values, which is finely 
gridded, at the specified modulation frequency. But the finer the gridding of the 
search region gets, the more correlations must be computed. Thus, a disadvantage 
with using a fine grid for searching over the hypothesized range of delay values is that 

20 the high cost of computing the correlations, especially when there is a considerable 
amount of sampled data, (i.e., the sampling duration is long) can be computationally 
prohibitive. Based on the foregoing, there is a clear need for a technique for estimating 
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the delay associated with a received signal in a cost effective fashion during signal 
processing. 

A very similar problem arises when the received signal is modulated by an un- 
known frequency and one needs to derive a fine-tuned estimate of this frequency from 
a set of IQ correlations obtained by searching over a coarse-grained set of frequencies 
representing the uncertain region. 
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SUMMARY OF THE INVENTION 

Techniques are provided for fine-tuning estimates of a delay value and other pa- 
rameters of interest associated with a sampled signal. One aspect of the invention is 
to use, for the sampled signal, calculations of the In Phase and Quadrature (I and 
5 Q) correlation integrals at a limited number of coarsely spaced points, wherein the 
calculations are performed for a set Ei of hypothesized delay values. These coarse- 
grained IQ correlations can be carried out by means of dedicated hardware, or the 
method described in U.S. Patent Application No. / , , filed on the same day 
herewith, entitled "Synthesizing Coherent Correlation Sums at One or Multiple Car- 
lo rier Frequencies Using Correlation Sums Calculated at a Coarse Set of Frequencies" , 
by inventors Anant Sahai and John Tsitsiklis, (Attorney Docket No. 60021-0012), or 
any other possible way. A subset of I and Q correlation values based on the coarse 
granularity calculations of the I and Q correlation functions is used to perform an 
interpolation to obtain fine-grained values of the I and Q integrals in the range of 
15 the delay, or other parameter, values of interest. Searching over such interpolated 
values effectively produces correspondingly fine precision estimates of the delay, or 
other parameter, values in the range of interest. 

Another approach used in this invention to calculate fine precision estimates is 
based on a template-matching approach whereby the received signal is compared 
20 with a parametric template. 

This invention can also be used to estimate parameters other than the delay, such 
as the modulation frequency or other characteristics of the received signal. 
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In general, these estimates can be given in the form of a single optimal value, an 
interval, a collection of intervals or a set of values. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention is illustrated by way of example, and not by way of limita- 
tion, in the figures of the accompanying drawings and in which like reference numerals 
refer to similar elements and in which: 
5 FIG. 1 is a block diagram that illustrates a system overview for processing a signal; 

FIG. 2 is a flowchart that illustrates a technique for fine-tuning estimates of delay 
in time associated with a received signal; 

FIG. 3 is a flowchart that illustrates a template-based approach to calculate fine- 
tuned estimates of the delay and other parameters; 
10 FIG. 4 is a block diagram illustrating a computer system on which embodiments 
of the invention may be implemented. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Techniques are provided for fine-tuning estimates of a delay and other parameter 
values associated with a sampled signal In the following description, for the purposes 
of explanation, numerous specific details are set forth in order to provide a thorough 
5 understanding of the present invention. It will be apparent, however, to one skilled in 
the art that the present invention may be practiced without these specific details. In 
other instances, well-known structures and devices are shown in block diagram form 
in order to avoid unnecessarily obscuring the present invention. 

SYSTEM OVERVIEW 

10 FIG. 1 is a block diagram that illustrates a system overview for processing a 
signal. System 100 comprises a signal source S, a receiver H, and a server D. In this 
example, the analog signal is transmitted to receiver H. The signal that is received 
at H is herein referred to as a "received signal". It is assumed that there is an 
unknown delay from the time the analog signal leaves signal source S and when it 

15 is received at receiver H. Such a delay is herein referred to as a "delay value" . In 
one embodiment of the invention, receiver H converts the analog signal into discrete 
sequence of values as a function of time by digitizing the received signal. The digitized 
received signal is herein referred to as a "sampled signal" or "sampled data". In 
one embodiment, receiver H transmits the sampled data to server D for processing. 

20 For the purpose of explanation, signal source S may be a transmitter or beacon 
like a Global Positioning System (GPS) satellite vehicle from a set of GPS satellite 
vehicles that is overhead the location of receiver H. Examples of receivers are Global 
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Positioning System receivers, cell phones with embedded signal receivers, Personal 
Digital Assistants (PDAs) with embedded signal receivers, etc. Further, in certain 
embodiments, receiver H may transmit the sampled data to a base station (not shown 
in FIG. 1), which in turn transmits the sampled data to server D for processing, as 
5 further described in U.S. Patent Application Serial No. 09/730,324, filed December 
4, 2000, entitled "Methods and Apparatus for Determining Location Using a Thin- 
Client", by inventors Benjamin Van Roy, Andrew Chou and Sherk Chung 

FUNCTIONAL AND OPERATIONAL OVERVIEW 
This invention reduces the amount of computation needed to calculate fine preci- 

10 sion estimates of the delay (or other parameters such as the exact Doppler shift) of 
the signal by means of a two-stage approach. Although the description here focuses 
on estimating delay, one of ordinary skill in the art should see how to use this inven- 
tion to estimate other parameters of interest. The first stage considers only a coarse 
set of possible delay values and produces one or more approximate delay estimates. 

15 The strength or likelihood of candidate delay values is evaluated either directly by the 
square-magnitude of the IQ correlations, I 2 + Q 2 , or by comparing the IQ correlations 
with a template (template-matching method). The IQ correlations are standard In- 
phase and Quadrature correlations of the sampled signal with a suitably modulated 
reference signal representing the known signal being searched for. 

20 The second stage refines these approximate delay estimates by performing either 
an interpolation or finer template-matching operation in a neighborhood of these 
coarse estimates. It avoids the expense of recalculating the IQ correlations at the 
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finer resolution. 

DETERMINING FINE-TUNED ESTIMATES OF DELAY VALUES 
FIG. 2 is a flowchart that illustrates the simplest embodiment of a technique 
for fine-tuning estimates of delay values associated with a received signal. FIG. 3 
5 illustrates the technique based on template-matching. In certain embodiments of the 
invention, server D performs the steps as described herein with reference to FIG. 2 
and FIG. 3. 

Coarse IQ Correlation: At block 220, coarse-grained calculations of IQ correla- 
tion integrals are performed for the sampled data at a set Ei of hypothesized delay 

10 values. The determination of a range of hypothesized delay values can be based on 
the uncertainty in certain basic variables like space and time, or might be known 
from other sources depending on the application at hand. Determining this range is 
further described in U.S. Patent Application Serial No. / , , filed on the same day 
herewith, entitled "Determining the Spatio-Temporal and Kinematic Parameters of a 

15 Signal Receiver and Its Clock by Information Fusion", by inventors Anant Sahai, An- 
drew Chou, Wallace Mann and Stefano Casadei, (Attorney Docket No. 60021-0014). 
In one embodiment of the invention, described further in U.S. Patent Application No. 
/ , , filed on the same day herewith, entitled "Signal Acquisition Using Data Bit 
Information", by inventors Anant Sahai, Wallace Mann, Andrew Chou and Benjamin 

20 Van Roy (Attorney Docket No. 60021-0011), the coarse-grained calculations are car- 
ried out together by dividing the sampled signal into a plurality of blocks of data. 
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Alternatively, the coarse-grained calculations can be carried out individually in the 

direct manner well known to those skilled in the art. 

The coarse-grained IQ correlations can also be carried out by means of dedicated 

hardware, or the approximate method described in U.S. Patent Application No. / , 
5 , filed on the same day herewith, entitled "Synthesizing Coherent Correlation Sums at 

One or Multiple Carrier Frequencies Using Correlation Sums Calculated at a Coarse 

Set of Frequencies", by inventors Anant Sahai and John Tsitsiklis, (Attorney Docket 

No. 60021-0012), or by any other known technique. 

The underlying IQ correlations at hypothesized delay r are denoted by the complex 
10 signal z(t) = J(r) + jQ(r). The coarse-grained IQ correlations can be represented 

by a discrete sequence z[k) of complex numbers where the k-th point in the sequence 

is z(t) evaluated at r = kA where A denotes the inter-sample time. 

Magnitudes: At block 230, the magnitudes \z[k]\ = y/P[k] + Q 2 [k] of the coarse- 
grained IQ correlation integrals are calculated. The magnitude of the continuous-time 
15 correlation function \z(r) | is assumed to be large (e.g. above a threshold) or maximum 
when r is at or near the true delay. 

Computation of coarse delay estimate(s) (or "initial estimate of delay val- 
ues"): At block 240, one or more coarse values of the delay, herein referred to as 
"initial estimates of delay values," are selected by using one of the following methods. 
20 The set of selected values is denoted by Yf v One method is to choose the single value 
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of kA that maximizes |;z(A;A)|: 



S' x = {argmax fcA |2(feA)|}. 



Another method is to select an interval of values around the maximum. One 
skilled in the art can see how statistical methods can be used to determine the size 
5 of this interval given a specified confidence level and suitable statistical assumptions 
relevant to the exact situation at hand. 

Another method is to choose all the values kA for which \z(kA)\ is above a 
threshold b: 

E[ = {kA: \z{kA)\ >b}. 

10 One skilled in the art can see how statistical methods can be used to determine 
this threshold given a specified confidence level and suitable statistical assumptions 
relevant to the exact situation at hand. 

Fine-grained values of interest (or "range of delay values of interest"): At 

block 250, a set S 2 = of delay values in the neighborhood of the initial estimate of 
15 delay value is selected to represent the new range of delay values of interest. For 
instance, these values of interest can be chosen to be: 



where jA is one of the selected coarse values of the delay, (also referred to as "selected 
initial estimate of delay value") and P is the desired precision. 




(1) 
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Fine-grained interpolation At block 260, fine-grained values for the I and Q cor- 
relation integrals are interpolated by using a subset of the coarse-grained calculations 
of the I and Q correlation integrals. The fine-grained values for the I and Q correla- 
tion integrals can then be searched to give rise to fine-tuned estimated delay values. 

5 Various interpolation techniques may be used for interpolating the fine-grained and Q 
correlation integrals. As one skilled in the art knows, there are a variety of well known 
interpolation techniques and which one is suitable may vary from implementation to 
implementation. One such technique is bandlimited interpolation. 

If the IQ correlation integral z(r) is presumed to have bandwidth 1/2A or less 

10 viewing delay r as the virtual time variable, then the fine-tuned value of the IQ 
correlation z(a 2 ), o 2 € S 2 is computed by interpolating the A-spaced samples z(kA) 
by means of a sine weighting function: 

( \ ST n.A^ 5 ^^- A:A)/A] 

The cutoff value Ri at which the sum is truncated can be determined according to 
15 the accuracy desired, as one skilled in the art can readily see. The larger the desired 
accuracy, the large R\ should be made. 

Finally, block 270 selects a set E' 2 of fine-tuned estimates of delay values by using 
one of the same set of possible methods that could be used in block 240. 

Alternatively, the interpolation above could be carried out using the real valued 
20 \z(t)\ rather than the complex z(r) if the approximation loss is acceptable. 

TEMPLATE-BASED APPROACH 
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FIG. 3 is a flowchart that illustrates another technique for fine-tuning estimates 
of delay values and other parameters associated with a received signal. For this 
technique, it is assumed that the ideal shape of the IQ correlation integral z(t) is 
known and is represented by a parametric template S^t) parametrized by /j,. When 
5 the hypothesized source produces and transmits a signal, the IQ correlation z[k] is 
assumed to contain an image of the template 5 M given by S^(t - a), where a is the 
delay. The observed IQ correlation integral z(t) is modeled to be 

z(t) = S^t -a) + noise. 

The parameters of the template (other than the delay), denoted collectively by 
10 the symbol account for the variability in the received signal due to a variety 
of modeled phenomena, which may include signal attenuation, frequency shifts, and 
multipath propagation. For example, the IQ correlation integral for GPS signals in the 
presence of multipath contains the interaction of two or more low-passed triangular 
waveforms that are scaled and shifted with respect to each other by amounts that vary 
15 in a continuum of possible values. Other variability is due to some signal processing 
parameters (such as the bandwdith of certain filters) that are known for the particular 
signal being processed but can vary between signals. 
One form for the template that is considered: 

S„{t -a) = A,Sf\t - a) + . . . + A N S { e N) (t - a), 

20 where the parameter /j, comprises both the linear parameters A±,...,A N and the 
remaining nonlinear parameters (other than the delay a) 9. 
Attorney Docket: 60021-0013 -13- 



The template-matching approach consists in finding estimates of the parameters 
for which the total weighted square error 

k 

is small or minimal. The weights iu„, which are optional and could just be set to 1, 
5 can be used to emphasize certain portions of the template which are known to be 
more reliable or important for the specific situation at hand. 

The description of the template-matching approach is as follows. 
At block 310 of FIG. 3, a parametric family of templates S^t), ji e M is chosen. 
If weights are used, a weighting function for each parameter fi e M is constructed. 
10 Block 310 can be executed off-line. More details about this block are in the section 
entitled "Template construction" . 

Block 320 determines a coarse set of delay values (also referred to as "initial 
estimates of delay values"), Ei, and non-linear parameters, 0i, which presumably 
contain a good approximation to the unknown value of the parameter 0 and delay 
15 a. A loop is initiated that hypothesizes each of these values and control is passed to 
Block 330. 

Block 330 performs a linear regression to estimate the linear parameters A = 
(Ax, ... , A N ) that minimize the weighted square error for each of the hypothesized 
non-linear parameter 0 X and delay values kA. For instance, if there is only one linear 
20 parameter, 

A(kA ft) = E,^(^A-iA)5 gl -(^A-iA)#1 ^ 
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Block 340 selects a subset of the hypothesized coarse values of the delay (also 
referred to as "selected initial estimates of delay values") and other nonlinear param- 
eters according to the value of the minimum total weighted square error A(fcA,0i), 
equal to the total weighed square error A M (fcA) minimized over all linear parameters 
5 A. 

One method is to choose the value of these parameters that minimizes A(fcA, 
Another method is to select a range of values around the minimim and represent 
them by an uncertainty interval. One skilled in the art can see how a statistical 
method can be used to determine the size of this interval given a specified confidence 
10 level and suitable statistical assumptions relevant to the exact situation at hand. 

Another method is to choose all the values for which A (A; A, #i) is below a threshold. 
One skilled in the art can see how a statistical method can be used to determine 
the threshold given a specified confidence level and suitable statistical assumptions 
relevant to the exact situation at hand. 
15 The selected subsets of coarse delays (also referred to as "selected initial estimates 
of delay values") and nonlinear parameters is denoted E' x and 0'^ 

Block 350 selects a fine-tuned range of delay values that are also referred to as 
a "range of delay values of interest" , and nonlinear parameters in the neighborhood 
of each selected coarse value (i.e., of each "selected initial estimate of delay values"). 
20 For example, the following fine-grained values of the delay (i.e., "range of delay values 
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of interest") can be considered: 

£ 2 = j^P + A;A, i = 0,±l,±p fcAeEi}. 

Block 360 is analogous to block 330, but is performed on the fine-grained set of 
delay values (range of delay values of interest) and nonlinear parameters. However, 
5 since the hypothesized values are much closer to each other, one can obtain significant 
computational savings by using an interpolation technique whenever the underlying 
functions are sufficiently smooth. 

For example, in the case when there is only one linear parameter A, and a ban- 
dlimited assumption is valid, then an interpolation method similar to the one used in 
10 block 260 can be used: 



fc 



If the bandlimited assumption is not valid, then one needs to solve the linear 
regression problem for each of the hypothesized fine-grained values and search over 
the results. For example, 

This computation is done in slightly different ways depending on whether the 
template // and the weights are encoded by an analytical expression or by a table. 
One skilled in the art can implement a system to evaluate the expression (2) and 
search for an estimate in a satisfactory way according to the relevant constraints and 
20 requirements by means of exhaustive search, divide-and-conquer, gradient descent, 
Newton's method, and a host of other optimization techniques. 
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One skilled in the art will also see that in a small neighborhood of a nominal value 
for a nonlinear parameter, the effect of that parameter on the template can be well 
approximated by a linear term using a Taylor expansion. Such a method can be used 
to help aid the search over some of the otherwise nonlinear parameters in small local 
5 neighborhoods. One skilled in the art can use this to speed up the search even further 
by linearizing nonlinear parameters based on the neighborhoods in question. 

For certain types of templates, it is possible to estimate the delay and amplitude 
of the template by means of an exact iV-point fitting of the function S^t - a) to the 
observed IQ correlation integral. This is accomplished by solving the system of N 
10 nonlinear equations 

S^kA - a) = z(kA), k = + N, 

whose unknowns are the delay a, and the N - 1 unknown parameters in \i of the 
template. One skilled in the art will notice that the although the above equations 
are written for a sequential set of correlation integrals, the TV equations can also be 
15 made from any set of N observed correlation integrals. 

An important application of this method is to searching for an estimate for the 
exact modulation frequency parameter J. The analogous calculations are performed, 
but with the search over the delay parameter a replaced with a search over the 
modulation frequency /. The equations are 

20 S^kA - /) = z(feA), k = f, f + 1, 

where A is now a frequency interval and the template S^f) can be taken equal to 
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the expression 



cm A sin2[7Tf/A] 



In this particular choice of template for estimating the frequency, the unknown param- 
eters are just the amplitude A and the / representing the true modulating frequency. 
5 This template is an approximation to the case where the underlying correlation sums 
are taken over a large number of data points. Other templates might be suitable if 
the underlying correlation sums are computed differently. 
Block 370 is exactly analogous to block 340. 

TEMPLATE CONSTRUCTION 
10 In one embodiment of this invention, which calculates delay estimates of a GPS 
signal, the basic template is obtained by convolving the appropriately lowpass filtered 
PRN code (i.e. the appropriate PRN code filtered by the lowpass filter H w>a ) with 
itself. This basic template will be denoted gw,a 

S^t -a) = Ag WA {t-a). 

15 The parameter W denotes the bandwidth of the filter and a denotes other parameters 

that describe the shape of the filter. The template has roughly the shape of a rounded 

triangular peak with a few smaller peaks on the sides. 

Among the causes that affect the shape and bandwidth of the filter H w ,a are 

the duration of the correlation integrals, the properties of the A/D converter (and 
20 filters therein), and the nature of the approximations during the calculations of the 

correlation integrals (As those in U.S. Patent Application No. / , , filed on 
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the same day herewith, entitled "Signal Acquisition Using Data Bit Information", 
by inventors Anant Sahai, Wallace Mann, Andrew Chou and Benjamin Van Roy, 
(Attorney Docket No. 60021-0011) and U.S. Patent Application No. / , , filed 
on the same day herewith, entitled "Synthesizing Coherent Correlation Sums at One 
5 or Multiple Carrier Frequencies Using Correlation Sums Calculated at a Coarse Set 
of Frequencies" , by inventors Anant Sahai and John Tsitsiklis, (Attorney Docket No. 
60021-0012)). 

In more sophisticated embodiments designed to achieve robustness to multipath 
phenomena, a composite template is constructed out of two or more basic templates. 
10 For example, a template with two components is obtained by summing one basic 
template corresponding to the direct path of propagation and one associated with a 
reflected path: 

S„(t -&) = A igw , a {t -a) + A 2 gwA t ~ ° ~ °')- 

Here, a' is the delay of the reflected waveform relative to the direct waveform and 
15 A u A 2 represent the complex amplitudes of the two components of the signal. 

In general, the parameters that characterize the templates are of two different 
types: specified parameters and searched parameters. For the two-composite tem- 
plates just described, W,a are specified parameters and Ai,A 2 ,<J,a' are searched 
parameters. The specified parameters are known or selected in advance of the appli- 
20 cation of this invention, whereas the searched parameters and the image parameters 
are estimated by means of the processing described by the system in FIG. 3. 
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Whereas the delay is often the main parameter of interest, many embodiments of 
the invention can utilize also one of more of the other estimated parameters. The 
exact Doppler shift of the signal is one other such parameter. One skilled in the art 
will see that with estimating Doppler, an added consideration is the duration of the 
5 received signal that is being correlated. The longer the duration, the more finely 
spaced in Doppler one needs to have even the coarse-grained points for the estimates 
to be reliable. One skilled in the art will see that the application of the template 
fitting approach is not limited to extracting parameters of interest. The goodness of 
the fit or likelihood can also be used to determine whether the modeled phenomena 
10 accurately explain the observed data. For example, even with no explicit modeling 
of multipath phenomena, the presence of multipath can sometimes be detected by 
noticing that the goodness of the fit is poor. In such cases, the overal system can take 
advantage of this information by discounting estimates that come from data which is 
known to be poorly modeled. 
15 In some embodiments of the invention, the weights io M are not used (or equivalently 
set equal to one). In others, the weights are chosen so as to emphasize the more 
reliable or significant portions of the template. One skilled in the art will notice that 
the early portion of a basic template is more reliable because it is less affected by 
reflected paths. In some embodiments, the weighting function w^t) is set equal to 
20 1 for t < t 0 , and equal to 0 for t > t 0 . The details of how w^t) is chosen can be 
determined according to the particular situation, as one skilled in the art can see. 
In some embodiments of the invention, an analytical expression of the signature 
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template and/or the weighting function is used. In certain cases, these analytical 
expressions can be used to fit directly the template to the received signal. Alter- 
natively, one skilled in the art can use these analytical expressions to implement an 
input/output algorithm that returns the value of the template and/or the weighting 

5 function at particular values. In other embodiments, the template and weighting 
functions are pre-computed and stored in a tabular database that can significantly 
reduce the amount of time required to carry out the calculations. One skilled in the 
art can easily determine whether an analytical or a tabular approach is more ap- 
propriate to the situation at hand, and such decision can even be made at run-time 

10 automatically by suitable automated logic. 

HARDWARE OVERVIEW 
An embodiment of the invention may be implemented using a computer system 
that includes a processor for processing information. The Computer system also in- 
cludes a main memory, such as a random access memory (RAM) or other dynamic 
15 storage device, coupled to a bus for storing information and instructions to be exe- 
cuted by the processor. The main memory also may be used for storing temporary 
variables or other intermediate information during execution of instructions to be ex- 
ecuted by the processor. The computer system further includes a read only memory 
(ROM) or other static storage device coupled to the bus for storing static information 
20 and instructions for the processor. A storage device, such as a magnetic disk or optical 
disk, is provided and coupled to the bus for storing information and instructions. 

The invention is related to the use of the computer system for implementing the 
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techniques described herein. According to one embodiment of the invention, those 
techniques are implemented by the computer system in response to the processor ex- 
ecuting one or more sequences of one or more instructions contained in main memory. 
Such instructions may be read into the main memory from another computer-readable 

5 medium, such as the storage device. Execution of the sequences of instructions con- 
tained in the main memory causes the processor to perform the process steps described 
herein. One or more processors in a multi-processing arrangement may also be em- 
ployed to execute the sequences of instructions contained in the main memory. In 
alternative embodiments, hard-wired circuitry may be used in place of or in combi- 

10 nation with software instructions to implement the invention. Thus, embodiments of 
the invention are not limited to any specific combination of hardware circuitry and 
software. 

The term "computer-readable medium" as used herein refers to any medium that 
participates in providing instructions to the processor for execution. Such a medium 

15 may take many forms, including but not limited to, non-volatile media, volatile me- 
dia, and transmission media. Non-volatile media includes, for example, optical or 
magnetic disks, such as the storage device. Volatile media includes dynamic memory, 
such as the main memory. Transmission media includes coaxial cables, copper wire 
and fiber optics, including the wires that comprise the bus. Transmission media can 

20 also take the form of acoustic or light waves, such as those generated during radio 
wave and infrared data communications. 

Common forms of computer-readable media include, for example, a floppy disk, 
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a flexible disk, hard disk, magnetic tape, or any other magnetic medium, a CD- 
ROM, any other optical medium, punch cards, paper tape, any other physical medium 
with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any 
other memory chip or cartridge, a carrier wave as described hereinafter, or any other 
5 medium from which a computer can read. 

Various forms of computer readable media may be involved in carrying one or more 
sequences of one or more instructions to the processor for execution. For example, 
the instructions may initially be carried on a magnetic disk of a remote computer. 
The remote computer can load the instructions into its dynamic memory and send 
10 the instructions over a telephone line using a modem. A modem local to computer 
system can receive the data on the telephone line and use an infrared transmitter to 
convert the data to an infrared signal An infrared detector coupled to the bus can 
receive the data carried in the infrared signal and place the data on the bus. The bus 
carries the data to the main memory, from which the processor retrieves and executes 
15 the instructions. The instructions received by the main memory may optionally be 
stored on the storage device either before or after execution by the processor. 

The computer system also includes a communication interface coupled to the bus. 
The communication interface provides a two-way data communication coupling to a 
network link that is connected to a local network. For example, the communication 
20 interface may be an integrated services digital network (ISDN) card or a modem to 
provide a data communication connection to a corresponding type of telephone line. 
As another example, the communication interface may be a local area network (LAN) 
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card to provide a data communication connection to a compatible LAN. Wireless links 
may also be implemented. In any such implementation, the communication interface 
sends and receives electrical, electromagnetic or optical signals that carry digital data 
streams representing various types of information. 

5 The network link typically provides data communication through one or more 
networks to other data devices. For example, the network link may provide a con- 
nection through the local network to a host computer or to data equipment operated 
by an Internet Service Provider (ISP). The ISP in turn provides data communication 
services through the worldwide packet data communication network now commonly 

10 referred to as the "Internet". The local network and the Internet both use electrical, 
electromagnetic or optical signals that carry digital data streams. The signals through 
the various networks and the signals on the network link and through the communi- 
cation interface, which carry the digital data to and from the computer system, are 
exemplary forms of carrier waves transporting the information. 

15 The computer system can send messages and receive data, including program code, 
through the network(s), the network link and the communication interface. In the 
Internet example, a server might transmit a requested code for an application program 
through the Internet, the ISP, the local network and the communication interface. 
In accordance with the invention, one such downloaded application implements the 

20 techniques described herein. 

The received code may be executed by the processor as it is received, and/or 
stored in the storage device, or other non-volatile storage for later execution. In this 
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manner, the computer system may obtain application code in the form of a carrier 
wave. 

In the foregoing specification, the invention has been described with reference to 
specific embodiments thereof. It will, however, be evident that various modifications 
5 and changes may be made thereto without departing from the broader spirit and scope 
of the invention. The specification and drawings are, accordingly, to be regarded in 
an illustrative rather than a restrictive sense. 
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